Alternatives to Eigenstate Thermalization 
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An isolated quantum many-body system in an initial pure state will come to thermal equilibrium if it satisfies 
the eigenstate thermalization hypothesis (ETH). We consider alternatives to ETH that have been proposed. We 
first show that von Neumann's quantum ergodic theorem relies on an assumption that is essentially equivalent to 
ETH. We also investigate whether, following a sudden quench, special classes of pure states can lead to thermal 
behavior in systems that do not obey ETH, namely, integrable systems. We find examples of this, but only for 
initial states that obeyed ETH before the quench. 
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Introduction and summary of ETH: Consider an isolated 
quantum iV-body system with Hamiltonian H. Let \a) de- 
note an eigenstate of H with eigenvalue E a , and let A de- 
note a few-body observable. The eigenstate thermalization 
hypothesis (ETH) states that (1) the diagonal matrix elements 
A aa = (a\A\a) change slowly with the state, with the differ- 
ence between neighboring values A a+1 Q+1 — A aa exponen- 
tially small in N, and (2) that the off-diagonal matrix elements 
A a p = (a\A\fi), a ^ (3, are themselves exponentially small 
in TV [1, 2]. ETH is suggested by various results in quan- 
tum chaos theory (in particular, Shnirelman's theorem [3] and 
Berry's random-wave conjecture [4]) for systems that have a 
chaotic classical limit. ETH has been verified numerically in 
a wide variety of quantum many-body systems that are suf- 
ficiently far (in parameter space) from points of integrability 
[5-8], but it certainly does not hold in systems that are inte- 
grable or near integrable [5-9] . 

Let |^>(t)) be the quantum state of the system at time r, 
given by the evolution of some initial state \tpj) = J2 a Ca \ce), 



-iHr/h 



e- iE " T/h \a) 



(1) 



with J2 a |C«| 2 = 1- The energy of the system is E = 
\C a \ 2 E a , and the quantum energy uncertainty is AE, 
where (AE) 2 = J2 a \C a \ 2 {E a - E) 2 . We assume that AE 
is algebraically small in N, (e.g., AE ~ N~ 1 / 2 E), as is the 
case for any state of a macroscopic system that could be real- 
istically prepared in a laboratory. (This is also true for sudden 
quenches between Hamiltonians with short-range interactions 
[5].) The time-dependent expectation value of A is 



(i(r)) = (V(r)|i^(r)} 



(2) 



a/3, 



and the long-time average of (A(t)) is 

A = lim - [ dt(A(T))=Y u \C a \ 2 A aa , (3) 

in the absence of degeneracies (which are not expected to oc- 
cur in chaotic systems without extra symmetries [10]). The 



right-hand side of Eq. (3) effectively sums \C a \ 2 A aa over an 
energy window of width AE that is centered on E. Accord- 
ing to ETH, A aa is approximately constant over this window. 
Thus, up to algebraically small corrections, the right-hand side 
of Eq. (3) has the same value as the microcanonical average 
of A over the same window, and is independent of the detailed 
pattern of values taken by the \C a \ 2 coefficients. Thus ETH 
results in the equality of time averages and thermal averages 
for a very broad class of initial states. 

From ETH, each term in the second sum in Eq. (2) is ex- 
ponentially small in N. However, the number of terms in this 
sum is exponentially large, and so, if the phases line up co- 
herently, the second sum can rival the first. If this happens 
at one particular time (say, r = 0), it will fail at sufficiently 
large later times, as the time-dependent phases go out of align- 
ment. This dephasing mechanism accounts for the approach 
to thermal equilibrium for an initially out-of-equilibrium state 
[2, 5]. 

The quantum ergodic theorem: In 1929, von Neumann 
proved a mathematical result which has been dubbed the 
quantum ergodic theorem (QET) [11]. An exegesis of it has 
been given by Goldstein et al. (hereafter GLTZ) [12]. GLTZ 
summarize QET, or "normal typicallity" as it has been more 
recently known, as follows: "for a typical finite family of com- 
muting macroscopic observables, every initial wave function 
from a microcanonical energy shell so evolves that for most 
times r in the long run, the joint probability distribution of 
these observables obtained from \ip(r)) is close to their mi- 
crocanonical distribution" [12]. More specifically, QET states 
that (A(t)) will be close to the microcanonical average of A 
in the following sense: \(A(t)) - (A) mc \ 2 < e 2 (A 2 ) mc forall 
but a fraction S of times t, where the subscript mc denotes the 
microcanonical average over the energy window of all states 
with nonzero C a , and e and 5 are small numbers. The proof 
of the theorem requires that all energy eigenvalue differences 
E a — Ep be nondegenerate, and an additional condition that 
is deemed technical by GLTZ, their Eq. (17). Here we point 
out that this condition is equivalent to ETH. Hence, von Neu- 
mann's proof of QET relies on ETH [13]. 

We follow the exposition of GLTZ, and consider the system 
to have a Hilbert space H with an exponentially large but fi- 
nite dimension D. We focus on a single observable A. GLTZ 
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partition the Hilbert space into "macrospaces" H u , with di- 
mension d v . Each H v is spanned by the eigenstates of A 
with eigenvalue a in a particular range centered on a value 
a v . Let P v be the projection operator onto % v , and consider 
its energy-basis matrix elements (a\P u \/3). The condition that 
must be assumed to prove QET is that, for each v, the off- 
diagonal elements (a ^ (3) must be exponentially small, and 
the diagonal elements (a\P u \a) must be exponentially close 
to f v = d v /D, the fraction of states in H v [11, 12]. 

We now argue that this condition is effectively equivalent 
to ETH. Consider the operator A cg = ~^2 U a v P v . This is 
a coarse-grained version of A itself, which can be written 
as A ~ ^2 a a\a)(a\, A cg is the same, but with the indi- 
vidual values of a replaced by their average values in each 
macrospace. The intention of von Neumann and GLTZ is 
that A cg should be indistinguishable from A in practical ex- 
periments. Next, consider the energy-basis expectation value 
(a | A C g | a). Since (a|Pi/|o;) is (by the GLTZ technical con- 
dition) exponentially close to f v , then from the definition of 
A cg we get (a|j4 cg |a) = ^„ a v f u , up to exponentially small 
corrections. Up to the additional small errors introduced by 
the coarse graining (which are assumed to be negligible), this 
last expression is equal to the trace of A, which in this simpli- 
fied model is to be identified with the microcanonical average 
of A. This is equivalent to the ETH statement that an energy- 
basis expectation value is equivalent to a microcanonical aver- 
age. Similarly, if (a|P„|/3) is exponentially small for 
so is (a|^4|/3), which is the other key part of ETH [2, 5]. 

Another way to phrase the equivalence is to note that both 
the GLTZ technical condition and ETH rely on exponential 
smallness of the overlap between an energy eigenstate and an 
eigenstate of an observable A that exhibits thermal behavior. 
As already noted, ETH can be justified by various results from 
quantum chaos theory. Thus, ETH provides a physical basis 
for the technical condition needed by QET. This results in a 
unification of two formerly disparate schools of thought on 
the foundations of statistical mechanics. 

Quantum quenches: As mentioned in the introduction, 
ETH has been shown to be satisfied in a variety of noninte- 
grable quantum systems. It has been found to break down 
only as one approaches integrable points [6-8], or in special 
regimes that are dominated by finite size effects, e.g., close 
to the atomic limit [14, 15]. Thermalization itself has been 
shown to be robust in nonintegrable systems after a (sudden) 
quench, once again, failing to occur close to integrable points 
[6, 16, 17] or the atomic limit [15, 18], and in localized disor- 
dered systems [19]. Here, by (sudden) quench, we mean that 
the system is prepared in an eigenstate of some initial Hamil- 
tonian (not necessarily the ground state) and then at r = the 
Hamiltonian is changed. 

We now consider systems that do not satisfy ETH; in partic- 
ular, we consider systems for which A aa varies significantly 
[that is, by an amount that is 0(N )] with a. It is easy to 
construct such systems; for example, any set of noninteract- 
ing degrees of freedom, with A corresponding to any one- or 
few-body observable, is in this class. Interacting systems that 



are integrable (with as many conserved charges as degrees of 
freedom) are in this class [20]. Equation (3) still applies to 
such systems; but now whether or not A is close to the ther- 
mal average of A depends strongly on the initial state, which 
is specified by the C a coefficients. If the values of the co- 
efficients \C a \ 2 in Eq. (3) provide an unbiased sampling of 
the matrix elements A aa , then we can expect A to be alge- 
braically close to the microcanonical average of A over the 
energy window specified by the AE of the initial state. 

The above scenario, however, does not occur in quenches 
between integrable systems, i.e., when the initial state is an 
eigenstate of an integrable system, and the time evolution 
is studied after changing some parameters in the Hamilto- 
nian while keeping the system integrable. Studies of several 
models have shown that A remains different from the ther- 
mal expectation as one approaches the thermodynamic limit 
[9, 21-34]. Even some special initial states that were seen to 
lead to A similar to the ones predicted in thermal equilibrium 
[9, 21, 28], have been recently shown not to result in the ther- 
malization A in the thermodynamic limit [32, 34]. 

Here we identify a class of initial states that leads to ther- 
mal behavior after a quench to an integrable point. The initial 
states we consider are eigenstates of an initial Hamiltonian 
Hj that is nonintegrable. This Hamiltonian is constructed by 
breaking the integrability of the final Hamiltonian Hp; i.e., we 
set Hj = Hp + "integrability breaking terms." The idea here 
is that the integrability breaking terms in the initial Hamilto- 
nian generate eigenstates that are unbiased combinations of 
eigenstates of the integrable (final) Hamiltonian. This is what 
leads to chaotic behavior as one departs from an integrable 
point [7], and ultimately allows ETH to be valid in noninte- 
grable systems. Hence, such initial states enable the desired 
unbiased sampling that does not occur in quenches between 
integrable systems. 

In order to show that this is indeed the case, we have stud- 
ied one-dimensional lattice systems of hard-core bosons and 
spinless fermions with the Hamiltonian 



L< 



V 



(hj - i) (h j+1 - |) + V' (hj - i) (h j+2 - |) 



, (4) 



where cj (£j) stands for the creation (annihilation) operator 

for hard-core bosons and fermions, hj = CjCj is the site occu- 
pation operator, L is the number of lattice sites, and t (t 1 ) and 
V (V) are the nearest (next-nearest) neighbor hopping and in- 
teraction, respectively. The t' , V' terms are the ones that make 
this Hamiltonian nonintegrable. We consider periodic bound- 
ary conditions, and the full diagonalization of the Hamiltonian 
is done using its translational symmetry. The filling is always 
taken to be N = L/3. 

Our quench protocol is then as follows, we generate an ini- 
tial state that is an eigenstate of Eq. (4) with t = V = 1 
(this sets our energy scale), t' — V ^ 0, and which lies 
within the sector of zero total momentum. The final Hamil- 
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FIG. 1. Time evolution of 8nt, after a sudden quench, for hard-core 
bosons (left panels) and spinless fermions (right panels) in lattices 
with L — 21 and L — 24 (N = L/S), and T = 3. Results are 
presented for quenches from t' = V' = 0.04 (a),(e), t' = V" = 0.16 
(b),(f), = V = 0.32 (c),(g), and t' = V = 0.64 (d),(h), to 
t' = V' = 0. In all cases t = V = 1 before and after the quench. 



tonian, after the quench, still has t = V = 1, but we 
set t' = V' = 0. This Hamiltonian is integrable. We 
then study the relaxation dynamics of various observables, 
as well as their description after relaxation, for many dif- 
ferent initial states. These initial states are eigenstates of 
Hamiltonians with different values of t' , V, and are selected 
so that the system can have different final effective temper- 
atures. The effective temperature T is calculated using the 
standard procedure for the canonical ensemble, i.e., such that 
E = Z- 1 lv{H F e-^"/ kBT ), where E = (ip^HFltpi) is the 
energy of the time evolving state, Z = Tr(e~ HF / kBT ) is the 
partition function, and fee (set to one in what follows) is the 
Boltzman constant. 

We have studied four observables, the kinetic energy K, the 
interaction energy U, the momentum distribution function hf. 
(which is the Fourier transform of the one-particle correlations 
pij = c\cj), and the structure factor Nk (which the Fourier 
transform of the density-density correlations Nij = fiihj). K 
and U are local observables while hk and Nk are nonlocal, 
and K and hk are single-body observables while U and Nk 
are two-body observables. In all cases studied, we found a 
similar qualitative behavior in those four quantities. Hence, 
we will only report results for hk, which is the one with the 
closest connection to ultracold gases experiments [35-38]. 

We are first interested in understanding how, after the 
quench, observables relax (if they do) to the long-time aver- 
age in Eq. (3). For hk, this can be conveniently quantified by 
calculating the normalized integrated difference 

Snk{T) = sjy))-^. 



(5) 



Results for Srik (r) are shown in Fig. 1 for bosons (left col- 
umn) and fermions (right column) and for two different sys- 
tem sizes. For very small quenches [Fig. 1(a), (e)], the initial 
state has a nonzero overlap only with very few close-by eigen- 
states of the final Hamiltonian. Because of this, the long-time 
average is close to the initial nk, and large oscillations (rel- 
ative to the time average) can be seen for both bosons and 




FIG. 2. Relative difference between as predicted by the micro- 
canonical ensemble and the long-time average, as a function of t' , V' 
in the initial Hamiltonian. Results are reported for hard-core bosons 
(a) and spinless fermions (b), for lattices with 24 sites (main panels) 
and 21 sites (insets). In all systems, we considered five different ini- 
tial states that resulted in five effective temperatures after the quench 
(T = 2, 3, 5, 7, and 10 in the plots). 



fermions. By increasing the amplitude of the quench, the 
long-time average becomes increasingly different from the 
initial momentum distribution, and the time fluctuations de- 
crease. This is because more eigenstates of the final Hamil- 
tonian are involved in the dynamics and dephasing becomes 
more efficient. It is interesting to note the strikingly large 
finite size effects seen for the spinless fermions [Fig. 1(e)- 
(g)], where the dynamics changes (improves) dramatically by 
increasing the system size from L — 21 to L ~ 24 [6, 7]. 
Overall, one can conclude from those results that, as the sys- 
tem size increases, nk relaxes to the long-time average rather 
quickly (r ~ h/t), and the time fluctuations around that aver- 
age become very small. Similar results were found for other 
observables and effective temperatures. 

We then are left to check how accurate statistical ensem- 
bles are when predicting TCk~. For our small finite systems, we 
use the microcanonical ensemble, and calculate the following 
normalized integrated difference to quantify its accuracy 



Efc l(^fc)mc ~n k \ 

Ek^ 



(6) 



The width AE of the energy window in the microcanonical 
ensemble is taken such that the results are robust to small 
changes of AE (in our systems AE ^0.1-0.3). 

Figure 2 depicts our results for Ank, calculated for bosons 
(a) and fermions (b), and for L = 24 (main panels) and 
L = 21 (insets). In that figure, it is apparent that nk is 



4 




FIG. 3. Normalized IPR (see text) for the same quenches depicted in 
Fig. 2. Results for hard-core bosons (fermions) are shown in the left 
(right) panels, and for L = 21 (L = 24) in the top (bottom) panels. 



nonzero, and large, for small values of t',V'. This is true 
even in the absence of a quench, when t' = V' = in the ini- 
tial and final Hamiltonian, and reflects the failure of ETH due 
to integrability. However, as the value of t' ', V is increased, 
rife is seen to decrease in all cases. This occurs for bosons and 
fermions at all effective temperatures and system sizes. For 
all observables and effective temperatures that we have stud- 
ied, we have found that Arife decreases with increasing system 
size. This, together with the understanding of the role of the 
integrability breaking terms in the initial Hamiltonian, sup- 
ports our claim that initial states that satisfy ETH (eigenstates 
of a nonintegrable Hamiltonian) will lead to thermalization in 
integrable systems, despite the fact that the latter do not satisfy 
ETH. 

In Fig. 2, one can also see that, in many instances, An./; 
for a fixed system size reaches a minimum value and then in- 
creases as t' and V' are increased. This occurs because, for 
large values of t', V', the initial state starts having large over- 
laps with eigenstates outside the microcanonical energy win- 
dow. Hence, even though the initial state is sampling more 
eigenstates of the final Hamiltonian, the energy density be- 
comes broad and the microcanonical ensemble once again be- 
comes a bad approximation for the long-time average. The 
fact that more eigenstates of the final Hamiltonian are part of 
the initial state can be quantified by means of the inverse par- 
ticipation ratio (IPR) IPR = 1/ J2 a \Ca\ 4 - This quantity has 
been shown to increase in eigenstates of nonintegrable many- 
body Hamiltonians as one departs from an integrable point 
[7]. On the other hand, the fact that the weight of the ini- 
tial state within the microcanonical window decreases if the 
value of t', V becomes too large can be quantified calculat- 
ing W — J2 a , \C a '\ 2 , where only eigenstates a' inside the 
microcanonical window are added. 

Then, one can compute a "normalized IPR," which is the 
product of the IPR and W. For finite systems, this quantity 
can tell us how effective t', V are in sampling states within 
the microcanonical window. Results for this quantity are de- 
picted in Fig. 3, for the same quenches depicted in Fig. 2. 



Figure 3 shows that, in the region where Arife exhibits a sharp 
decrease in Fig. 2, the normalized IPR increases. In addition, 
where Arife saturates or increases in Fig. 2, the normalized 
IPR saturates or decreases in Fig. 3. This allows one to un- 
derstand the overall behavior of Arife in Fig. 2. However, we 
should stress that the fact that t' ', V (or whatever other term is 
used to break integrability) cannot be made too large is only 
a concern for finite systems. As long as these terms are kept 
O(N ) and the interactions have finite range, the energy width 
of any initial state after the quench will vanish in the thermo- 
dynamic limit [5], and the initial state will only sample states 
within the microcanonical window. 

Conclusions: In this Letter we have further probed the 
role of the eigenstate thermalization hypothesis (ETH) as the 
key dynamical feature of systems that come to thermal equi- 
librium. We have shown that von Neumann's quantum er- 
godic theorem (QET) relies on a technical assumption that 
is in fact essentially equivalent to ETH. We have also exam- 
ined whether, after a sudden quench, eigenstates of the initial 
Hamiltonian can lead to thermal behavior in systems that do 
not obey ETH. We have found that this is possible, but only 
if those eigenstates obeyed ETH before the quench. These re- 
sults further support the fundamental role of ETH in thermal 
behavior of quantum many-body systems. 
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cussion on QET. 



[1] J. M. Deutsch, Phys. Rev. A 43, 2046 (1991). 

[2] M. Srednicki, Phys. Rev. E 50, 888 (1994); J. Phys. A 29, L75 

(1996); J. Phys. A 32, 1163 (1999). 
[3] A. I. Shnirelman, Usp. Mat. Nauk 29, 181 (1974). 
[4] M. V. Berry, J. Phys. A 10, 2083 (1977). 
[5] M. Rigol, V. Dunjko, and M. Olshanii, Nature 452, 854 (2008). 
[6] M. Rigol, Phys. Rev. Lett. 103, 100403 (2009); Phys. Rev. A 

80, 053607 (2009). 
[7] L. F. Santos and M. Rigol, Phys. Rev. E 81, 036206 (2010); 

Phys. Rev. E 82, 031130 (2010). 
[8] C. Neuenhahn and F. Marquardt, arXiv: 1007.5306 (2010). 
[9] A. C. Cassidy, C. W. Clark, and M. Rigol, Phys. Rev. Lett. 106, 

140405 (2011). 

[10] O. Bohigas, M.-J. Giannoni, and C. Schmit, Phys. Rev. Lett. 52, 
1 (1984). 

[11] J. von Neumann, Zeitschrift fur Physik 57, 30 (1929), translated 
to English in European Phys. J. H 35, 201 (2010). 

[12] S. Goldstein, J. L. Lebowitz, R. Tumulka, and N. Zangh, Euro- 
pean Phys. J. H 35, 173 (2010). 

[13] von Neumann actually proved that this "technical" condition 
holds for "most" macroscopic observables. However, "most" 
was defined by a probability distribution that treats each ob- 
servable as a random matrix in the energy basis. For real world 
observables, which are physically motivated and not explicitly 
built as random matrices, that condition has not been proved. 

[14] G. Roux, Phys. Rev. A 81, 053604 (2010). 

[15] M. Rigol and L. F. Santos, Phys. Rev. A 82, 01 1604(R) (2010). 



5 



[16] S. R. Manmana, S. Wessel, R. M. Noack, and A. Muramatsu, 

Phys. Rev. Lett. 98, 210405 (2007). 
[17] M. C. Bafiuls, J. I. Cirac, and M. B. Hastings, Phys. Rev. Lett. 

106, 050405 (2011). 
[18] C. Kollath, A. M. Lauchli, and E. Altaian, Phys. Rev. Lett. 98, 

180601 (2007); G. Biroli, C. Kollath, and A. M. Lauchli, Phys. 

Rev. Lett. 105, 250401 (2010). 
[19] C. Gogolin, M. P. Muller, and J. Eisert, Phys. Rev. Lett. 106, 

040401 (2011). 

[20] B. Sutherland, Beautiful Models (World Scientific, Singapore, 
2004). 

[21] M. Rigol, V. Dunjko, V. Yurovsky, and M. Olshanii, Phys. Rev. 
Lett. 98, 050405 (2007); M. Rigol, A. Muramatsu, and M. Ol- 
shanii, Phys. Rev. A 74, 053616 (2006). 

[22] M. A. Cazalilla, Phys. Rev. Lett. 97, 156403 (2006). 

[23] P. Calabrese and J. Cardy, J. Stat. Mech. p. P06008 (2007). 

[24] M. Cramer, C. M. Dawson, J. Eisert, and T. J. Osbome, Phys. 
Rev. Lett. 100, 030602 (2008). 

[25] T. Barthel and U. Schollwock, Phys. Rev. Lett. 100, 100601 
(2008). 

[26] M. Eckstein and M. Kollar, Phys. Rev. Lett. 100, 120404 
(2008). 

[27] M. Kollar and M. Eckstein, Phys. Rev. A 78, 013626 (2008). 



[28] D. Rossini, A. Silva, G. Mussardo, and G. E. Santoro, Phys. 
Rev. Lett. 102, 127204 (2009); D. Rossini, S. Suzuki, G. Mus- 
sardo, G. E. Santoro, and A. Silva, Phys. Rev. B 82, 144302 
(2010). 

[29] A. Iucci and M. A. Cazalilla, Phys. Rev. A 80, 063619 (2009); 

New J. Phys. 12, 055019 (2010). 
[30] D. Fioretto and G. Mussardo, New J. Phys. 12, 055015 (2010). 
[31] J. Mossel and J.-S. Caux, New J. Phys. 12, 055028 (2010). 
[32] P. Calabrese, F. H. L. Essler, and M. Fagotti, Phys. Rev. Lett. 

106, 227203 (2011). 
[33] M. A. Cazalilla, A. Iucci, and M.-C. Chung, arXiv: 1106.5206 

(2012). 

[34] M. Rigol and M. Fitzpatrick, Phys. Rev. A 84, 033640 (2011). 
[35] T. Kinoshita, T. Wenger, and D. S. Weiss, Nature 440, 900 
(2006). 

[36] S. Hofferberth, I. Lesanovsky, B. Fischer, T. Schumm, and 

J. Schmiedmayer, Nature 449, 324 (2007). 
[37] C.-L. Hung, X. Zhang, N. Gemelke, and C. Chin, Phys. Rev. 

Lett. 104, 160403 (2010). 
[38] S. Trotzky, Y.-A. Chen, A. Flesch, I. P. McCulloch, 

U. Schollwock, J. Eisert, and I. Bloch, arXiv:1101.2659 (2011). 



